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ABSTRACT 


The  effect  of  increased  atmosphere  rotation  rate  on  low  perigee 
altitude  orbits  and  the  effect  of  underestimation  of  atmosphere  rotation 
rate  on  orbit  determination  are  displayed  using  simulation  results.  Cross  - 
track  changes  in  satellite  position  due  to  increased  rotation  rate  are  small. 
Intrack  changes,  though  larger,  are  small-error  sources  in  orbit  determina¬ 
tion  and  short-term  prediction  if  a  seven-parameter  fit  is  used.  The  angles 
between  satellite  inertial  velocity  and  wind  vectors  that  produce  zero 
tangential  acceleration  and  maximum  normal  acceleration  are  derived  for 
any  wind  and  for  the  special  case  of  a  circular  orbit  in  a  rotating  atmosphere. 
This  analysis  explains  the  increased  decay  rates  of  some  near-polar  orbits 
due  to  a  rotating  atmosphere  and  the  inability  to  predict  this  effect  with 
frequently  used  approximations. 
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SECTION  I 


INTRODUCTION  AND  SUMMARY 


The  King-Hele  (Reference  1)  results,  displaying  high  values  of  atmos¬ 
phere  rotation  rate  (1.0  to  1.  6  times  earth  rate),  prompted  a  simulation 
study  to  determine  the  effect  of  a  large  rotation  rate  (twice  earth  rate)  on 
the  prediction  of  satellite  position  for  low  perigee  altitude  orbits.  Initially 
it  was  believed  that  large  crosstrack  errors  in  polar  orbit  determination 
and  short-term  prediction  could  result  from  assuming  an  atmosphere  rotation 
rate  smaller  than  the  true  rate. 

The  Aerospace  Corporation  orbit  determination  program,  TRACE,  was 
used  in  the  two-part  simulation  study.  In  the  first  part,  radial,  intrack,  and 
crosstrack  differences  were  obtained  by  comparing  pairs  of  orbits.  Initial 
conditions  were  identical  for  each  orbit  of  a  pair,  but  one  was  subjected  to 
an  atmosphere  rotating  at  twice  earth  rate,  the  other  to  an  atmosphere 
rotating  at  earth  rate.  Crosstrack  position  changes  due  to  a  doubled  rate 
proved  to  be  small,  even  for  polar  orbits.  However,  unexpected  intrack 
changes  (increased  decay  rates)  occurred  for  polar  orbits. 

The  second  part  of  the  study  showed  that  intrack  perturbations  due  to 
a  doubled  atmosphere  rotation  rate  are  largely  predictable  if  a  seven-parameter 
fit  is  used  in  orbit  determination  with  the  ballistic  coefficient  W/(C^A)  as  the 
seventh  parameter.  For  each  orbit,  simulated  observation  data  were  obtained 
for  the  orbit  in  an  atmosphere  rotating  at  twice  earth  rate.  A  seven-parameter 

,  f  t  ,  ,  :  .  t  i  r  t  r  f  t  t  t  l  f  '  '  ' 

fit  to  twelve  hours  of  simulated  tracking  data/W&s  then  obtained.  For  this  fit 
and  for  the  subsequent  prediction,  the  atmosphere  waft„assumed  to  rotate  at 
earth  rate.  Using  the  resulting  parameter  values,  position  was  then  predicted 
up  to  twelve  hours  beyond  the  fit  span.  Radial,  intrack,  and  crosstrack  dif¬ 
ferences  were  obtained  by  comparing  the  predicted  positions  with  those  of  the 
true  orbit.  All  resulting  differences  were  small. 
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In  the  simulations,  substantial  increases  in  the  decay  rates  of  polar 
orbits  were  not  expected  from  commonly  used  approximations.  To  explain 
this,  the  angle  between  satellite  inertial  velocity  and  wind  vectors  that 
produces  zero  tangential  acceleration  due  to  wind  is  derived.  It  is  slightly 
less  than  90  deg;  smaller  angles  decrease  decay  rates  and  larger  angles 
increase  them.  The  angle  for  maximum  normal  acceleration  is  also  derived; 
it  is  slightly  larger  than  90  deg.  Equivalent  results  are  obtained  for  the 
special  case  of  circular  orbits  in  a  rotating  atmosphere.  The  equations  yield 
the  inclination  for  zero  intrack  acceleration  and  for  maximum  crosstrack 
acceleration  due  to  atmosphere  rotation  for  circular  orbits.  The  inability  to 
predict  these  results  with  the  approximate  equations  in  References  2  and  3 
is  discussed  for  the  circular-orbit,  rotating-atmosphere  case. 
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SECTION  II 


SIMULATION  RESULTS 


A.  EFFECTS  OF  ATMOSPHERE  ROTATION 
RATE  ON  ORBITS 


Initial  interest  was  in  the  effects  of  increased  atmosphere  rotation 
rate  on  polar  orbits.  When  some  effects  proved  to  be  substantial,  interest 
turned  to  the  implications  with  respect  to  orbit  determination  and  short¬ 
term  prediction  techniques;  these  are  discussed  In  "Section  II.  B.  This 
section  presents  results  for  equatorial  orbits  as  well,  to  indicate  the 
variation  of  effects  with  orbit  inclination. 

For  each  of  the  cases  studied,  the  satellite  position  at  15-minute 
intervals  was  obtained  for  two  orbits  over  a  12-hour  period.  Initial  con¬ 
ditions  were  identical  for  both  orbits.  Force  models  used  were  also  identical 
but  with  one  exception:  For  one  orbit  the  atmosphere  rotated  at  earth  rate; 
for  the  other  it  rotated  at  twice  earth  rate.  Subtracting  appropriate  satellite 
position  coordinates  of  the  former'from  those  of  the  latter  at  corresponding 
times  yielded  the  radial,  intrack,  and  crosstrack  differences.  Keplerian 
element  differences  were  also  obtained. 

Table  I  shows  the  extremes  of  the  departures  of  the  differences  from 
zero  after  12  hours  (approximately  8  revolutions).  These  extremes  occurred 
during  the  eighth  revolution  in  all  cases.  The  minimum  and  maximum  cross - 
track  differences  simply  indicate  the  occurrence  of  substantial  extremes  in 
both  positive  and  negative  directions  (due  primarily  to  a  monotonic  inclination 
change).  Other  differences  grew  monotonically  except  for  minor  variations 
during  each  revolution. 

The  quantities  shown  are  perigee  altitude  h  ,  apogee  altitude  h  , 

P  a 

inclination  i,  latitude  of  perigee  6^,  and  semimajor  axis  a.  The  ARDC  1959 
atmosphere  and  a  spherical  earth  were  used;  and  W/(C^A)  was  taken  as 
80  lb  /ft2. 
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Table  I.  Extreme  Differences  after  Eight  Revolutions 
Due  to  Doubling  Atmosphere  Rotation  Rate 


Case 

Orbit 

Differences 

h  (n  mi) 
P 

h  (n  mi) 

SL 

i  (deg) 

6p  (deg) 

Radial 

(ft) 

Intrack 

(ft) 

Crosstrack 

(ft) 

a 

(ft) 

(deg) 

1 

80  > 

90 

-60 

1179 

min,  -165 

-31 

-0.00045 

180 

34 

max,  1 66 

.2 

80 

90 

-70 

1416 

min,  -186 

-39 

-0.00056 

180 

0 

max,  186 

3 

70 

90 

-168 

3089 

min,  -391 

-81 

-0.0011 

230 

0 

max,  393 

4 

80 

0 

1864 

-39932 

o' 

1092  ! 

0 

180 

0 

5 

70 

0 

3857 

-73065 

0 

1924 

0 

230 

0 

Cases  1  to  3  involve  polar  orbits.  Cases  4  and  5  involve  equatorial 
orbits  which,  except  for  inclination,  are  like  those  of  Cases  2  and  3,  respec¬ 
tively.  These  pairs  of  cases  display  the  variation  of  effects  with  inclination. 
Case  1,  with  perigee  latitude  at  34  degrees  but  otherwise  the  same  as  the 
polar  orbit  of  Case  2,  shows  the  effect  of  reduced  winds  at  perigee  due  to 
increased  perigee  latitude. 

The  possibility  of  large  crosstrack  differences  was  the  primary  concern 
initially;  however,  it  is  clear  from  Cases  1  to  3  that  crosstrack  differences 
are  minor.  (Maximum  crosstrack  acceleration  occurs  for  inclinations  slightly 
greater  than  90  deg.  See  Section  IV.  C. )  It  is  easily  shown  that  the  polar 
orbit  crosstrack  differences  of  Table  I,  which  take  place  near  the  poles,  very 
nearly  correspond  to  the  inclination  changes  in  the  last  column;  the  perturbed 
orbit  is  to  tue  right  of  the  comparison  orbit  at  one  pole,  to  the  left  at  the 
other  pole. 
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The  relatively  large  intrack  differences  for  polar  orbits  were  unexpected. 
It  is  not  obvious  intuitively  that  increased  wind  normal  to  the  orbit  plane  will 
produce  the  increased  decay  rate  indicated  by  the  radial,  intrack,  and  semi¬ 
major  axis  differences  for  Cases  1  to  3.  By  comparing  Case  2  with  Case  4 
and  Case  3  with  Case  5,  the  polar  orbit  intrack  difference  is  seen  to  be  about 
4  percent  of  the  equatorial  orbit  intrack  difference  and  oppositely  directed. 

(The  equatorial  orbit  difference  is  the  maximum  possible  difference  for 
posigrade  orbits,  since  the  rotating  atmosphere  wind  is  greatest  at  the 
equator  and  the  wind  vector  lies  in  the  orbit  plane. )  That  is,  there  is  a 
significant  loss  of  satellite  energy  due  to  increased  wind  normal  to  the  polar 
orbit  plane. 

It  is  shown  in  Section  IV.  B  that  for  circular  orbits  (horizontal  velocity 
vectors),  the  inclination  for  zero  intrack  acceleration  due  to  atmosphere 
rotation  is  slightly  less  than  90  deg.  (Increased  decay  rates  occur  for  larger 
inclinations,  decreased  decay  rates  for  smaller  inclinations.)  From  Eq.  (37), 
this  inclination  is  computed  to  be  approximately  88.  3  deg  at  the  equator  for  an 
orbit  like  that  of  Case  3  with  an  atmosphere  rotating  at  twice  earth  rate.  The 
Case  3  90-deg  inclination  thus  differs  by  only  1.  7  deg  at  perigee  from  that 
which  would  produce  no  change  in  decay  rate  due  to  twice  the  earth  rate  atmos¬ 
phere  rotation,  indicating  sensitivity  of  the  phenomenon  to  inclination  change. 
(Commonly  used  approximations  in  the  literature  do  not  predict  these  effects. 
See  Section  V. ) 

B.  EFFECTS  OF  UNDERESTIMATED  ATMOSPHERE  ROTATION 

RATE  ON  ORBIT  DETERMINATION  AND  PREDICTION 

Tracking  was  simulated  for  a  12-hour  period  for  three  orbits  in  an 
atmosphere  rotating  at  twice  earth  rate.  Then,  by  assuming  an  atmosphere 
rotating  at  earth  rate,  seven  parameter  fits  of  the  tracking  data  were  obtained 
with  the  ballistic  coefficient  W/(C^A)  as  the  seventh  parameter.  The  result¬ 
ing  epoch  parameter  values  were  used  to  generate  predicted  positions  up  to 
12  hours  past  the  12-hour  tracking  interval.  True  satellite  coordinates  were 
then  subtracted  from  corresponding  predicted  coordinates  to  obtain  radial, 
intrack,  and  crosstrack  prediction  errors. 
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Results  are  displayed  in  Tables  II,  ID,  and  IV.  Epoch  parameter 

values  shown  in  Table  II  are  for  right  ascension  a,  geocentric  latitude  6, 

angle  between  velocity  vector  and  geocentric  vertical  (3,  azimuth  of  velocity 

vector  from  true  north  Az,  distance  from  geocenter  r,  magnitude  of  velocity 

vector  V  ,  and  the  ballistic  coefficient  W/(C_.A). 
o  u 

The  three  orbits  are  similar  in  shape  with  perigee  altitude  at  about 
68  n  mi  and  apogee  altitude  at  about  209  n  mi.  Orbit  1  is  a  polar  orbit  with 
perigee  at  34°N  latitude.  Orbits  2  and  3  have  inclinations  of  65  and  30  deg, 
respectively,  with  perigee  at  the  equator  (where  maximum  wind  occurs). 

In  the  simulations,  the  ARDC  1959  atmosphere  and  the  standard  TRACE 
geopotential  model  were  used.  Tracking  data  consisted  of  exact  range, 
azimuth,  and  elevation  at  30 -second  intervals  during  periods  of  visibility  for 
six  tracking  stations.  The  numbers  of  sets  of  range,  azimuth,  and  elevation 
-were  69,  49,  and  57  for  Orbits  1,  2  and  3,  respectively.  In  orbit  determina¬ 
tion,  the  standard  deviations  used  were  100  feet  for  range  and  0.  1  degree 
for  azimuth  and  elevation. 

Table  II  shows  true  epoch  parameter  values  for  the  three  orbits  and 
the  "corrections"  required  in  orbit  determination  due  to  assuming  an 
atmosphere  rotating  at  earth  rate.  Note  the  substantial  correction  in 
W/(CjpA)  for  Orbit  3,  which  experiences  the  largest  component  of  wind  in 
the  direction  of  satellite  motion  because  of  its  low  inclination.  This  correction 
indicates  that  in  the  orbit  determination  inplane  forward  acceleration  (due  to 
unmodeled  high  wind)  is  attributed  to  reduced  drag.  The  largest  of  the  absolute 
values  of  observation  residuals  for  the  converged  seven-parameter  fits  are 
shown  in  Table  III. 

The  prediction  capability  of  the  epoch  parameter  values  for  these  fitted 
orbits  is  indicated  in  Table  IV;  it  shows  the  minimum  and  maximum  prediction 
errors.  Errors  were  obtained  at  15-minute  intervals  from  epoch  plus 
12  hours  to  epoch  plus  24  hours. 

The  orbits  of  this  section  are  not  the  same  as  those  of  Table  I,  in  which 
however,  it  is  clear  that  the  intrack  component  of  satellite  position  can  change 
bymanynr'les  over  a  24-hour  period  if  the  atmosphere  rotation  rate  is  doubled. 
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Table  III.  Maxima  of  Absolute  Values  of  'uals  Using  Exact  Observations 
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The  Table  IV  results  indicate  that  if  a  seven-parameter  orbit  determination 
is  used  and  an  atmosphere  rotating  at  earth  rate  is  assumed,  even  a  doubled 
atmosphere  rotation  rate  does  not  lead  to  substantial  short-term  prediction 
errors. 

C.  SIMULATION  CONCLUSIONS 


•  Crosstrack  errors  in  orbit  determination  and  short-term 
prediction  due  to  underestimation  of  atmosphere  rotation 
rate  are  small,  even  for  the  exaggerated  underestimation 
used  in  the  simulations. 

•  Intrack  effects  on  most  low-altitude  orbits  (including  90  deg 
inclination  orbits)  due  to  increased  atmosphere  rotation  rate 
are  substantial.  However,  if  a  seven-parameter  fit  is  used, 
these  effects  are  masked  and  produce  only  small  intrack 
errors  in  orbit  determination  and  short-term  prediction. 


SECTION  in 


EFFECTS  OF  A  GENERALIZED  WIND 


The  analysis  in  this  section  applies  to  any  wind.  It  is  used  in  the 
Section  IV  analysis  of  intrack  and  crosstrack  effects  of  a  rotating  atmosphere 
on  circular  orbits,  which  helps  explain  the  preceding  simulation  results. 

The  variation  of  tangential  acceleration  A^,  and  normal  acceleration  A^  due 
to  wind  is  discussed.  An  equation  is  derived  for  determining  whether  wind 
instantaneously  increases  or  decreases  satellite  velocity.  The  angle 
between  satellite  velocity  and  wind  vectors  for  maximum  normal  accelera¬ 
tion  is  obtained. 

Figure  1  shows  earth-centered  inertial  velocities  of  the  wind  W  and 
the  satellite  V  ,  satellite  velocity  relative  to  the  air  V_,  and  angles  between 
these  vectors,  y  and  0.  Accelerations  due  to  (the  no-wind  case)  and  V 
are  directed  opposite  to  these  velocities  and  have  the  magnitudes 


A  =  KV2 


(2) 


where  K  is  considered  the  same  for  both  velocities  in  this  derivation;  that  is, 
the  drag  coefficient  is  the  same  for  all  directions  of  satellite  motion  relative 
to  the  atmosphere. 

A.  TANGENTIAL  ACCELERATION 

A^,  is  investigated  first.  It  is  given  by 

A,p  =  A  cos  0  =  KV2cos  0  .  (3) 
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Figure  1.  Velocity,  Acceleration,  and  Wind  V 


From  Figure  1 ,  it  is  evident  that 


=  V  cos  0  =  V  -  W  cos  v 
o  1 

(4) 

=  V2  +  w2  -  2V  W  cos  y 

0  o  * 

(5) 

From  Eqs.  (3)  and  (4), 

At  =  KVVt  =  KV(Vo  -  W  cos  y)  •  (6) 


By  using  V  from  Eq.  (5),  this  may  be  written  as 

At  =  K(Vq  -  W  cos  y)  (V*  +  W2  -  2VqW  cos  y)1  /2  .  (7) 

The  behavior  of  AT  with  increasing  y  is  readily  obtained  from  Eq.  (7). 
The  terms  containing  cos  y  have  negative  signs.  It  follows,  for  fixed  Vq 
and  W,  that  A^,  increases  monotonically  as  y  goes  from  0  to  180  deg.  Further, 
it  is  clear  from  Figure  1  that  if  y  =  0  deg 

AT  =  A  <  Aq  (8) 

implying  increased  velocity  due  to  wind,  and  if  y  =  180  deg 

AT  =  A  >  Aq  (9) 


implying  decreased  velocity. 
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The  zero  of  A^,  -  Aq  is  now  investigated.  Introducing  VQ  -  W  cos  y 
into  the  radical  in  Eq.  (7)  and  multiplying  yields 


At  =  K  [v4  +  V2W2  -  (4V2W  +  2V0W3)  cos  y  +  (5V2W2  +  W4)  cos2y 

'}in  ' 


-  2V  W3  cos3y 
o 


A^,  is  the  rate  of  change  of  VQ <.  For  no  change  due  to  wind, 


(H) 


Now  we  let  yQ  be  the  value  of  y  for  which  Eq.  (11)  holds.  It  is  clear  that 

4 

Eq.  (11)  occurs  when  the  sum  of  the  terms  after  VQ  in  Eq.  (10)  is  zero.  The 
result  is  the  equation  in  cos  yQ 

2V  W2cos3v  -  (5VZW  +  W3)  cos2y  +  (4V3  +  2V  W2)  cos  y  -  V2W  =  0 
O  •  o  O  '  1  o  O  o  '  1  o  o 

(12) 

From  Eqs.  (8)  and  (9)  and  the  discussion  following  Eq.  (7),  it  follows  that 
only  one  value  of  y  between  0  and  180  deg  makes  A^,  =  Aq.  Further,  the 
alternation  of  signs  for  the  terms  of  Eq.  (12)  implies  that  cos  yQ  2  0.  Con¬ 
sequently,  0  £  yQ  £  90  deg.  For  example,  if  Vq  =  25,  000  ft/sec  and 
W  =  500  ft/sec,  yQ  =  89.  6  deg;  that  is,  if  y  <  89-  6  deg,  A^,  <  Aq,  implying 
increased  velocity  due  to  wind;  if  y  >  89.  6  deg,  A^  >  Aq,  implying  decreased 
velocity. 

Because  W  is  small  compared  with  Vo»  yQ  is  always  near  90  deg,  so 

that  cos  y  -  0.  Also,  for  any  realistic  value  of  W,  the  absolute  value  of  the 

°  3 

coefficient  of  cos  y  in  Eq.  (12)  is  much  larger  than  are  those  of  cos  y  and 
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2  3 

cos  y  because  of  the  term  4Vq.  A  good  approximation  for  yq  can  therefore 
be  obtained  from 

(4V3  +  2V  W2)  cos  v  a  V2W  .  (13) 

'  o  o  '  To  o  '  ’ 

2  3 

This  is  done  by  neglecting  2VqW  ,  which  is  small  compared  with  4Vo»  and  by 
using 


cos  \0*=  (nY2  -  y  )  sir/2  -  yQ 


(14) 


This  yields 


V 


o 


W 

4Vo 


(15) 


which  is  accurate  to  within  two  minutes  of  arc  for  the  example  in  the 
preceding  paragraph. 

It  follows  from  both  Eqs.  (12)  and  (15)  that  yQ  £  90  deg,  and  from  Eq.  (15) 
thatyo  decreases  nearly  linearly  with  increasing  W  for  realistic  values 
of  W.  Again,  if  y  <  YQ>  satellite  velocity  increases  compared  with  no  wind; 
if  y  >  yQ,  satellite  velocity  decreases.  Note  that  yQ  is  independent  of  K  and 
hence  of  atmospheric  density  and  ballistic  coefficient.  Equations  (12)  and  (15) 
provide  a  criterion  for  the  instantaneous  sign  of  the  change  in  satellite 
velocity  due  to  wind,  whatever  the  model. 

B.  NORMAL  ACCELERATION 

An  is  given  by 

An  =  KV2sin  0  =  KVVn  =  KVW  sin  y  .  (16) 
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If  Eq.  (5)  is  used  in  Eq.  (16), 


AN  =  KW  sin  y  (V^  +  W2  -  2VqW  cos  v)1 /2  .  (17) 

As  y  goes  from  0  to  180  deg,  the  radical  V  increases  monotouically.  If  the 
properties  of  sin  y  are  used,  it  is  then  clear  that  A^  varies  from  zero  at 
Y  =  0  deg  to  a  maximum,  and  then  to  zero  at  y  =  180  deg. 

Let  Ym  he  the  value  of  y  for  which  is  a  maximum.  By  rewriting 
Eq.  (17)  we  obtain 

AN  =  KW^(V2  +  W2  -  2VqW  cos  y)  (1  -  cos2y)J  1/2  .  (18) 

The  maximum  occurs  when  the  derivative  of  the  radicand  is  zero.  After 
simplification,  this  condition  is  expressed  as 

3V  W  cos2y_  -  (V2  +  W2)  cos  y  -  V  W  =  0  .  (19) 

o  m  o  m  o 

The  only  meaningful  solution  is 


By  inspection  it  is  seen  that  cos  Ym  -  0,  so  that 


ym  2  90  deg 


(21) 


For  example,  with  VQ  =  25,000  ft/sec  and  W  =  500  ft/sec,  Ym  =  91.2  deg. 
Note  that  Ym»  like  yq»  is  independent  of  atmospheric  density  and  ballistic 
coefficient. 


SECTION  IV 


ROTATING  ATMOSPHERE  EFFECTS 


Consider  a  circular  orbit  about  a  spherical  earth  with  a  spherically 
symmetric  atmosphere  rotating  at  a  rate  u  about  the  earth's  axis.  In  the 
following,  the  magnitudes  of  the  intrack  and  crosstrack  components  of 
velocity  relative  to  the  atmosphere  and  of  acceleration  due  to  the  atmosphere 
are  investigated.  Implications  concerning  noncircular  orbits  are  discussed. 

A.  VELOCITY  AND  ACCELERATION* VARIATIONS  '  - 

Figure  2  shows  inclination  i  and,  at  P,  the  satellite  latitude  6,  and  the 
angle  between  satellite  velocity  ard  wind  vectors  y.  In  the  spherical  right 
triangle  FQR, 


•in(f  -  \)  =  co.  y  =  fff-j  .  (22) 

Figure  3  shows  the  vectors  occurring  at  P  in  Figure  2.  Figure  3  is 
the  same  as  Figure  l,  except  that  tangential  and  normal  components  are  in 
this  case  intrack  (subscript  I)  and  crosstrack  (subscript  C)  components.  If 
r  is  the  satellite's  distance  from  the  geocenter, 

W  =  ru  cos  6  .  (23) 

a  .  . 

If  Eqs.  (22)  and  (23)  are  used  in  Eq.  (4), 

VT  =  V  cos  0  =  V  -  rw  cos  i  .  (24) 

I  o  a  17 


Also,  from  Figure  3, 


V £  =  V  sin  0  =  W  sin  y 


(25) 
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By  substituting  W  from  Eq.  (23)  and  sin  y  from  Eq.  (22)  in  Eq.  (25)  and 
simplifying, 

V_  =  V  sin  0  =  rw  (cos26  -  cos2i)*^2  .  (26) 

Note  from  Eq.  (24)  that  Vj  is  the  same  at  all  points  in  the  orbit  for  a  given 

inclination;  that  is,  (for  0  deg  <  i  <  90  deg),  the  tendency  to  decrease  Vj, 
produced  by  reduction  of  y  as  the  satellite  moves,  say,  northward  from  the 
equator,  is  exactly  countered  by  the  reduction  in  wind  with  increasing 
latitude.  For  retrograde  orbits,  the  tendency  to  increase  Vj  produced  by 
■"increase  of  y  as  the  satellite  moves  northward  is  also  countered  by  the 
reduction  in  wind.  In  contrast,  from  Eq.  (26),  has  a  maximum  at  the 
equator  (6  =  0,  y  =  i)  and  is  zero  at  maximum  latitude  ( |  6  I  =  i  or  180  deg  -  i, 
y  =  0  or  180  deg  ). 

Intrack  and  crosstrack  acceleration  magnitudes  are 

Aj  =  KVZcos  0  (27) 

Ac  =  KYZsin  0  .  (28) 


Although  Vj  is  constant  throughout  the  orbit,  Aj  is  only  approximately  con¬ 
stant.  This  can  be  seen  by  using  Eq.  (24)  in  Eq.  (27)  to  obtain 


KVj 

”  cos  0 


K 


cos  0 


(V  -  rw  cos  i) 
'  o  a 


(29) 


Thus  Aj  varies  inversely  as  cos  0  for  any  specified  inclination.  But  cos  0  31 

for  realistic  W  and  V  so  that  the  variation  of  A.  is  small.  In  fact, 

o  I 


At  as  K(V 
I  o 


rw  cos  i) 

SL 


(30) 
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is  a  fair  approximation.  Note,  however,  that  for  i  =  90  deg,  Eq.  (30)  is 
2 

Aj  ^  KV^,  indicating  no  wind  effect;  the  intrack  differences  for  Cases  1,  2, 
and  3  in  Table  I  are  an  indication  of  the  error  in  Eq.  (30). 

The  variations  of  Aj  and  A^,  with  latitude  and  inclination  are  now  con¬ 
sidered.  Using  Eqs.  (24)  and  (26)  in  Eqs.  (27)  and  (28)  yields 

AT  =  KVV.  =  KV(V  -  ru  cos  i)  (31) 

X  X  O  3 

Ac  =  KVVC  =  KVr«a(cos26  -  cos^)1  ^  .  (32) 

By  using  Eqs.  (22)  and  (23),  Eq.  (5)  can  be  written  as 

V  =  (v^  +  r2u2cos26  -  2V  ru  cos  i)*  ^2  .  (33) 

0  3.  0  3 

For  a  specified  inclination,  only  cos  6  varies  in  Eq.  (33),  so  that  V  is  a 
maximum  at  the  equator  and  a  minimum  when  the  satellite  is  at  maximum 
latitude  (north  or  south).  (This  can  also  be  concluded  from  the  variations  of 
Vj  and  discussed  previously.  )  If  a  similar  argument  is  used  for  the 
radicand  in  Eq.  (32),  it  follows  from  Eqs.  (31),  (32),  and  (33)  that 

(a)  For  any  inclination,  Aj  decreases  from  a  maximum  at  the 
equator  to  a  minimum  at  the  greatest  distance  from  the 
equator  (|6|  =  i  or  180  deg  -  i).  Again,  from  Eq.  (29),  the 
variation  is  seen  to  be  small. 

(b)  For  any  inclination,  Ar  decreases  from  a  maximum  at  the 
equator  to  zero  at  the  greatest  distance  from  the  equator. 

Next  consider  the  effect  on  Aj  and  A^.  of  variation  of  inclination,  holding 

6  fixed.  For  6*0,  it  is  clear  that  the  possible  range  of  values  of  i  is 

I  6 1  <  i  £  180  deg  -  |  6|,  For  any  such  subset  of  0  £  i  £  180  deg,  the  terms 

involving  cos  i  in  Eqs.  (31)  and  (33)  increase  monotonically  as  i  increases. 

2 

The  term  -cos  i  in  Eq.  (32)  increases  to  zero  at  i  =  90  deg,  then  decreases 
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thereafter  as  i  increases.  Again,  from  Eqs.  (31),  (32),  and  (33),  it 
follows  that 

(c)  For  any  latitude,  Aj  increases  monotonically  as  i  increases 
through  its  possible1  range. 

(d)  For  any  latitude,  Ac  increases  from  zero  (at  i  =  |6|)  to  a 
maximum  and  then  decreases  to  zero  (at  i  =  180  deg  -  |6|)  as 
i  increases  through  its  possible  range. 

B.  ZERO  INTRACK  ACCELERATION 

The  preceding  conclusions  (c)  and  (d)  are  closely  related  to  those  of 

Section  III.  First,  it  is  apparent  from  Figures  2  and  3  that  Aj  <  Aq  for 

i  =  0  deg  and  that  A^  >  Aq  for  i  =  180  deg.  These  conditions  indicate, 

respectively,  for  the  circular  orbit  and  rotating  atmosphere,  decreased 

decay  rate  and  increased  decay  rate  due  to  wind.  The  condition  A^  =  Aq 

(no  change  in  satellite  speed  due  to  wind)  was  investigated  in  Section  III.  A. 

The  approximation  Eq.  (13)  is  applied  here  by  using  Eqs.  (22)  and  (23)  with 

2 

iQ  as  the  inclination  that  gives  Aj.  =  Aq.  Again  if  2VqW  is  neglected, 
this  yields 


ru  cos  6 
a 


(34) 


(36) 


Substituting  in  Eq.  (35)  gives 


r u  cos 
.  it  a 
i  —  —  - 
1o"  2 


4V 


(37) 


Thus,  i  -  90  dee  and  i  depends  on  5;  that  is,  i  increases  as  latitude 
o  o  o 

increases.  Note  that  there  are  circular  orbits  with  inclinations  slightly 
below  90  deg  that  experience  increased  decay  rate  at  the  equator  and 
decreased  decay  rate  near  the  poles. 


C.  MAXIMUM  CROSSTRACK  ACCELERATION 


Simililarly,  if  Eqs.  (22)  and  (23)  are  used  in  Eq.  (20),  im,  the  incli¬ 
nation  for  maximum  A^,  then  is  given  by 


cos  ‘m  = 


V  +  * 
o  T  V 


7  '  (v» +  14wZ  +  £ ) 


1/2 


.  (38) 


Clearly  cos  i _  SO  implying  i _  >  90  deg.  Like  i  ,  i  depends  on  6,  since 

m  m  o  m  ^ 

W  =  ru  cos  6.  Since  i  as  y  when  y  as  90  deg,  and  if  the  results  of  the  example 

2k 

following  Eq.  (21)  are  used,  it  is  clear  that  im  is  close  to  90  deg.  Further, 
inspection  of  Eq.  (38)  reveals  that  i  approaches  90  deg  as  I  6  I  approaches 
90  deg,  since  W  approaches  zero. 

D.  COMPARISON  WITH  SIMULATION  RESULTS 

The  qualitative  statements  that  follow  relate  the  preceding  analytical 
results  to  the  results  in  Section  II.  A.  These  statements  are  made  with  the 
following  observations  in  mind:  First,  the  differences  in  Table  I  are  in  the 
same  direction  as,  but  smaller  in  magnitude  than,  those  that  would  have 
been  obtained  by  comparing  twice  earth  rotation  rate  and  no  atmosphere 
rotation  situations.  Second,  the  simulation  results  are  for  low -eccentricity 
orbits;  however,  radial  velocity  components  are  small,  particularly  in  the 
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vicinity  of  perigee  where  most  of  the  atmocphere  rotation  effects  are 
experienced.  The  simulated  effects  are  taken  as  resembling  the  circular 
orbit  effects  of  the  analysis. 

•  The  crosstrack  differences  for  the  polar  orbits  of  Cases  1  to 

3  in  Table  I  are  slightly  less  than  maximum;  maximum  cross- 
track  differences  would  occur  for  orbital  inclinations  slightly 
greater  than  90  deg  (see  Section  IV.  C). 

•  As  indicated  in  Section  II.  A,  the  90  deg  inclinations  are  at 
most  only  slightly  greater  than  the  instantaneous  inclination 
at  any  point  in  each  orbit  that  would  produce  zero  intrack 
acceleration.  This  is  shown  in  Section  IV.  B.  Again,  the 
intrack  differences  obtained  for  these  90-deg.inclination  orbits 
show  the  sensitivity  to  inclination  change  of  intrack  acceleration 
due  to  atmosphere  rotation. 

•  Intrack  differences  for  the  polar  orbits  of  Cases  1  to  3  in 

Table  I  are  an  indication  of  the  error  in  the  approximate  Eq.  (30), 
which  predicts  no  intrack  acceleration  for  polar  orbits  due  to 
atmosphere  rotation.  Eq.  (29)  does  predict  this  effect.  This  is 
discussed  further  in  Section  V. 

•  From  paragraph  (a)  following  Eq.  (33),  it  is  to  be  expected  that 
a  change  in  argument  of  perigee  will  produce  a  change  in  intrack 
acceleration  due  to  atmosphere  rotation.  This  is  illustrated  by 
the  intrack  differences  for  Cases  1  and  2  in  Table  I. 


SECTION  V 


NOTES  ON  THE  LITERATURE 

An  increased  decay  rate  of  polar  orbits  due  to  a  rotating  atmosphere 
is  not  expected  from  frequently  used  approximations  in  the  literature.  This 
led  to  the  analyses  in  Sections  III  and  IV.  The  inclination  separating  decreased 
and  increased  decay  rate  situations  is  given  approximately  by  Eq.  (37). 
Although  this  inclination  is  close  to  90  deg  (and  less  than  90  deg  except  at 
the  poles),  the  intrack  effect  on  a  90  deg  orbit  is  substantial  (see  Table  I). 
Further,  for  some  orbits  with  inclinations  slightly  less  than  90  deg,  a  rotating 
atmosphere  increases  the  decay  rate.  Thus,  any  generalization  implying  that 
all  posigrade  orbits  have  reduced  decay  rates  due  to  atmosphere  rotation  is 
misleading.  Discussion  here  of  two  references  explains  the  inability  of 
their  approximate  equations  to  yield  these  conclusions. 

In  Reference  2,  an  equation  is  used  that  can  be  written  for  a  circular 
orbit  as 


r  co  2  2 

V  =  V  -  ru>  cos  i  +  --- a  (sin  i  -  sin  6)  +  smaller  terms  (39) 
o  a  2V 

o 

using  the  notation  of  Section  IV.  All  terms  after  -r«  cos  i  are  then  neglected. 

a  2 

(Note  that  as  i  approaches  90  deg  the  term  involving  sin  i  becomes  larger 
than  that  involving  cos  i. )  The  resulting  approximation  is 


V  ~  V  -  ru  cos  i 
o  a 


When  applied  to  a  circular  orbit,  Eq.  (2.  10)  in  Reference  3  can  be  rewritten  as 


,,2  2  (.  rwa  ,\2  2  ; 

V  =  V°\  "  \T"  COS  lJ  +  r  “ 


2/  2*  2.. 

I  (cos  0  -  cos  1)  . 

cl 


2  2  2  2 

Reference  3  then  neglects  r  w  (cos  6  -  cos  i),  which  also  yields  Eq.  (40). 
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By  comparing  Eq.  (40)  with  (24),  the  approximations  in  both  references 
are  seen  to  be  expressions  for  Vjj  that  is,  the  right  member  of  Eq.  (40) 
requires  division  by  cos  6  to  make  the  equation  exact.  For  a  polar  orbit, 

Eq.  (40)  yields 


At  3  KV2cos  0 
I  o 

The  exact  value  from  Eq.  (29)  for  a  polar  orbit  is 

KVo 

AI  =  cos  0_ 


(42) 


(43) 


and  for  no  wind  the  exact  value  is 


Aj  =  KV2  .  (44) 

2 

Thus,  Eq.  (42)  differs  by  a  factor  of  cos'  0  from  the  exact  value  and  is 
even  smaller  than  the  exact  value  for  no  wind;  it  implies  a  decreased  rather 
than  an  increased  orbit  decay  rate  due  to  wind  for  a  polar  orbit. 

Comparing  Eqs.  (42),  (43),  and  (44)  leads  to  an  estimate  of  the  intrack 
error  that  results  if  Eq.  (42)  is  used  instead  of  Eq.  (43)  for  polar  orbits. 
Subtracting  Eq.  (44)  from  Eq.  (42)  yields  the  change  in  Aj  due  to  wind,  using 
the  approximation,  as 


KV.2cos  0  -  KV2  =  KV2(cos  0  -  1 ) 
o  o  o 


(45) 


Subtracting  Eq.  (44)  fr>im  Eq.  (43)  yields  the  exact  change  in  Aj.  due  to  wind,  as 


KV 


cos  0 


-  KV2  =  KV2  ( — -  l) 
o  o  Vcos  0  / 


(46) 
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•  ■muiri 


The  ratio  of  Eqs.  (45)  and  (46)  reduces  to  -cos  6,  which  is  approximately  -1. 
That  is,  for  polar  orbits,  the  acceleration  due  to  wind,  if  Eq.  (42)  is  used, 
has  about  the  same  magnitude  as  the  true  value,  but  is  oppositely  directed; 
or,  the  error  in  Eq.  (42)  is  almost  twice  the  exact  intrack  acceleration  due 
to  wind  for  polar  orbits. 

A  quantitative  implication  for  polar  orbits  is  obtained  if  Table  I  is 
considered  to  display  approximately  the  intrack  differences  between  orbits 
subjected  to  an  atmosphere  rotating  at  earth  rate  and  those  subjected  to  a 
nonrotating  atmosphere.  The  conclusion  of  the  preceding  paragraph  implies, 
to  first  order,  that  intrack  differences  using  Eq.  (42)  for  intrack  acceleration 
would  be  equal  to  the  intrack  differences  for  Cases  1,  2,  and  3  (polar  orbits) 
of  Table  I,  but  oppositely  directed. 

Further,  it  is  again  noted  that  the  Case  2  and  3  differences  are  approxi¬ 
mately  4  percent  of  the  Case  4  and  5  differences,  respectively.  It  follows 
that  the  errors  in  predicting  intrack  differences  due  to  wind  using  Eq.  (42) 
for  Cases  2  and  3  would  be  approximately  8  percent  of  the  intrack  differences 
due  to  wind  for  Cases  4  and  5.  Cases  4  and  5  are  selected  in  this  comparison 
because  these  equatorial  orbits  display  the  maximum  intrack  differences  due 
to  wind  for  all  similar  posigrade  orbits.  When  compared  with  these  maxi¬ 
mum  differences,  the  errors  of  Eqs.  (40)  and  (42)  in  predicting  intrack 

differences  due  to  wind  for  polar  orbits  are  seen  to  be  substantial. 

2 

The  error  in  V  ,  if  Eq.  (40)  is  used,  is  obtained  by  subtracting  values  of 

2 

V  from  Eqs.  (40)  and  (41),  the  latter  being  exact.  The  error  is 

A(V^)  =  r^w^(cos^5  -  cos^i)  .  (47) 

ct 

2 

That  is,  the  accuracy  of  intrack  acceleration  Aj  =  KV  cos  0  in  using 

Eq.  (40)  is  latitude-  and  inclination-dependent.  By  inspection,  it  is  seen 
2 

that  A(V  )  £  0;  the  true  relative  velocity  is  generally  greater  than  the  approxi- 

2 

mate  relative  velocity.  Also,  A(V  )  =  0  only  if  |6  |  =  i  or  180  deg  -  i.  The 

2  2  2 

maximum  value  of  A(V  ),  r  u>  ,  occurs  when  6  =  0  deg  and  i  =  90  deg.  Since 

ct  “ 
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the  variation  of  cos  6  is  small,  the  polar  orbit  intrack  acceleration  errors 
discussed  are  very  nearly  maximum.  Further,  the  error  does  not  change 
rapidly  as  i  departs  from  90  deg  (with  perigee  remaining  near  the  equator). 
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